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Ultracold neutrons (UCNs) can be bound by the potential of terrestrial gravity and a reflecting 
mirror. The wave function of the bound state has characteristic modulations. We carried out an 
experiment to observe the vertical distribution of the UCNs above such a mirror at Institut Laue- 
Langevin in 2011. The observed modulation is in good agreement with that prediction by quantum 
mechanics using the Wigner function. The spatial resolution of the detector system is estimated to 
be 0.7 /an. This is the first observation of gravitationally bound states of UCNs with submicron 
spatial resolution. 
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Terrestrial gravity is the most common force expe- 
rienced in everyday life. However, experimental mea- 
surements of quantum mechanical bound states in the 
Earth's gravitational field was only started in the last 
decade by the pioneering works of Nesvizhevsky et al. 
|1| using ultracold neutrons (UCNs). UCNs are neu- 
trons with kinetic energies lower than the Fermi pseudo- 
potential of suitable materials (~ 100 neV) and hence 
can be totally reflected by the material surfaces at any 
angle of incidence. The wave function ip(z) of a UCN 
in the terrestrial gravity field obeys the Schrbdingcr 
equation in the vertical direction z. The cigenstatcs 
of this system are linear combinations of Airy functions 
[2J]. The probability distribution in the vertical position 
of UCNs in the bound states, namely, the sum of the 
absolute squares of eigenfunctions, has a characteristic 
modulation. The eigenstate is specified by two scales, 
the length (h 2 / (2m 2 g)) = 5.87 fim and the energy 

(mg 2 ti 2 /2) = 0.602 peV where H is the reduced Planck 
constant, m is the neutron mass and g is the gravitational 
acceleration. There are therefore two ways to observe 
the bound state, to measure its energy or its position. 
Recently first measurements of the difference between 
the eigenenergies using the transitions of UCNs between 
quantum states in terrestrial gravity potential have been 
reported [3(- The ability of experiments to observe the 
spatial distribution of the bound state is limited by the 
spatial resolution (several microns) of current slow neu- 



tron detectors. 

To observe the spatial distribution of gravitationally 
bound states with high precision, wc developed a novel 
technique as shown in Figure Q] with its three main com- 
ponents [4|. The Cartesian coordinate system (x,y,z) is 
defined in Figure [T] (b). Z is the axis on the pixelatcd 
detector corresponding to the magnified height z. The 
collimating guide and magnification rod have a width of 
50 mm along y. As shown in Figure[T](b), incident UCNs 
pass through the collimating guide in which they settle 
into gravitationally bound states above the flat bottom 
mirror. The ceiling removes UCNs whose wave functions 
significantly penetrate the ceiling. The height distribu- 
tion of the surviving UCNs is magnified by a cylindri- 
cal rod which acts as a convex mirror. After reflection 
at the rod surface, UCNs are detected by a CCD-bascd 
pixelated detector. Using this setup, we performed an 
experiment to observe the spatial distribution of grav- 
itationally bound states during a period of 17 days in 
August 2011 at the Institut Laue-Langevin (ILL). De- 
tails of the experimental apparatus are discussed in the 
following. 

Since the statistics of the experiment is inevitably re- 
stricted by the narrow height of the guide, a high inten- 
sity UCN source is required. We used the UCN beam line 
PF2 @1 at ILL, the world's highest intensity steady UCN 
source. The horizontal velocity (v x ) distribution of UCNs 
was measured using standard time-of-flight technique. 
The measured distribution is nearly Gaussian with a 
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FIG. 1. Experimental setup, (a) general view and (b) the 
main components. The two thick bent arrows in (b) cor- 
respond to the trajectories of neutrons flying horizontally 
above the bottom mirror with the height of the bottom mirror 
z = fim and that of the ceiling z — 100 fim. The Roman 
numerals denote the calculation steps. 



mean of 9.4 m/s and standard deviation of 2.8 m/s. 

Inside the collimating guide made of glass with height 
h = 100 /im, the energy of UCNs is quantized. The 
ceiling of the guide removes UCNs with high vertical en- 
ergy due to its microscopic surface roughness. Once such 
a neutron is reflected by the ceiling, the large horizon- 
tal velocity component is converted into vertical velocity 
component. The chance to hit the ceiling again is then 
enormously enhanced. Numerous collisions cause UCN 
loss by absorption or upscattering. The arithmetic mean 
roughness is 0.4 fmx. In addition an absorptive mate- 
rial of Gd-Ti-Zr alloy (54/35/11) was deposited on the 
glass by vacuum evaporation |6| at the Kyoto University 
Research Reactor Institute (KURRI). The thickness of 
the layer is 200 nm and its potential is calculated to be 
-13.9 - 26.5i neV. 

At the end of the collimating guide, a cylindrical glass 
rod of 3 mm radius magnifies the distribution of UCN 
like a convex mirror. The geometrical arrangement of the 
magnification system is shown in Figured] (b). The distri- 
bution of 100 /im in height z is magnified to ~ 2.5 mm in 
the position Z on the detector, hence the average mag- 
nification power is about 25. The glancing angle was 
only 20° in order to make the critical energy of reflection 
high. Differences in v x cause dispersion of the parabolic 
trajectories and smear the distribution. This dispersion 
is estimated to be less than 0.1 /jm of the height z. The 



rod was precisely ground at the Crystal Optics Inc. and 
finely polished at the Research Center for Ultra-Precision 
Science and Technology, Osaka University. Furthermore 
a Ni layer of 200 nm was deposited on the polished glass 
surface, increasing the potential from 100 neV to 200 neV 
so that all UCNs exiting the guide were totally reflected. 
The Ni deposited surface has an arithmetic mean rough- 
ness of 1.9 nm, much smaller than the wavelength of UCN 
~ 100 nm. Hence, the diffused reflection from the surface 
roughness of the rod can be neglected. 

For high resolution two-dimensional detection, a 
charge-coupled device (CCD) was used. Since slow neu- 
trons would pass through the sensitive volume of CCDs 
without ionization, they must be converted to charged 
particles. To retain the intrinsic spatial resolution, a 
200 nm thin 10 B neutron converter was evaporated di- 
rectly on to a CCD. After neutron capture a and 7 Li 
particles are released. A back-thinned CCD |7| was the 
base of the detector. The pixel size is 24 /xm x 24 /im and 
the sensitive area is 24.576 mm x 6.000 mm. The length 
along the Z-axis is 6 mm. An incident charged particle 
creates electron-hole pairs inside the Si layer and loses 
energy. About 300,000 electrons per MeV are created 
and spread before reaching electrodes and being detected 
as a two-dimensional cluster. The barycenter of the de- 
posited charges corresponds to the incident position of 
the charged particle. The spatial resolution along the 
i?-axis was measured to be 3.35 ± 0.09 /im [8j]. 

The setup was installed inside a vacuum chamber to 
prevent neutrons from interacting with air. We evacuated 
the chamber to 10 Pa before the experimental run. The 
vacuum pump was disconnected during the measurement 
to reduce vibration. 

A neutron shutter Cd blades was installed inside an 
acrylic box as shown in Figure [T] (b) to shut off the UCN 
beam during the readout of the CCD. The box was con- 
nected to the beam pipe and the vacuum chamber by 
plastic bellows, to prevent the transmission of vibrations. 
The shutter box and the plastic bellows were filled with 
Helium gas to avoid scattering in air. 

The magnetic shield of mu-metal covering our experi- 
mental apparatus reduces the external magnetic field to 
about 1/100. As shown in Figure[T](a), the detector sys- 
tem was installed on granite and anti- vibration tables [9j 
to reduce vibration from the floor. 

The horizontality of the detector system was less than 
0.1 mrad and monitored by an inclinometer. The remain- 
ing effects due to external magnetic field, vibration and 
horizontality were estimated to be negligible. 

The distribution of UCNs on the pixelated detec- 
tor can be calculated using quantum mechanics. The 
state of a UCN can in general be written as a su- 
perposition of the n-th energy eigenstate as ^(z,t) = 
J2 n a n^n{z) exp(—iE n t / h) where a n satisfies £) n \a n \ 2 = 
1 and ip n (z) is normalized to give J dz tp n (z) = 1. The 
experimental result is an average over many incoherent 



UCNs. Since the phase of a n is randomly and uniformly 
distributed, the average of the absolute squared value of 



the superposition becomes |^(z)| 2 = J2 T 



l \^n{z)Y 



where the time dependence and interference terms are 
averaged out. In general, the quantum state of UCNs is 
treated as a mixed state. The state is described as a den- 
sity matrix, p = I]„Pn| - 0n)(V'n| where p n = \a n \ 2 is the 
probability of the n-th state and \tp n ) is the correspond- 
ing state vector. The calculation was in four steps (I, II, 
III, IV) from upstream toward the detector, as shown in 
Figure [T] (b). In the first three steps, the probability of 
the eigenstate p n at the end of the guide is calculated 
and then, in the last step, the UCN distribution on the 
pixelated detector is derived. Each step is discussed in 
the following paragraphs. 

(I) The state before the guide is a quantum state 
which can be approximated by a wave packet whose 
motion is described by classical mechanics because the 
de Broglie wavelength of 100 nm is much smaller than 
the geometrical scale. The energy distributions of the 
states before and after the transition at the guide en- 
trance is calculated from the allowed phase space area 
of the z direction in which UCNs can enter the guide 
between < z < 100 pm. The kinetic energy of UCN 
~ 100 ncV is much greater than the allowed vertical en- 
ergy of ~ 10 pcV, therefore the phase space density can 
be approximated to be uniform. Hence, the energy dis- 
tribution of the states at the guide entrance, namely, the 
probabilities of the states, is proportional to the allowed 
phase space area. It should be noted that the eigenstates 
with vertical energy larger than mgh (h = 100 pm) in- 
side the guide is different from that of without the ceiling, 
because the wave function becomes zero at the height of 
the ceiling. In the following, values with tildes denote 
that in the case with the ceiling at z = h. 

(II) For UCN loss inside the guide, we assume a phe- 
nomenological loss rate of the n-th state as T n + B n 
where T n and B„ arc respectively the loss by the ceil- 
ing and the bottom mirror. T n is assumed to be propor- 
tional to the probability of finding a UCN in the rough- 
ness region as r„ = 7 • f h _ 2 $dz \ipn(z)\ 2 where 7 is 
the constant for the loss and 5 = 0.4 pm is the arith- 
metic mean roughness of the ceiling (26 is the average 
width of the roughness region) [10J. B n is assumed to 
be proportional to the bouncing number per unit time 
of a bouncing motion above a floor in classical mechan- 



ics as B„ = P ■ (g /2\/2)\J m/ E n where (3 is the constant 
for the loss. Hence, the probability at the guide exit p n 
is written in terms of that at the guide entrance £» ra (0) 
as p n oc Pn(O) (exp (—l/v x ■ (r n + B„))) v where I is the 
length of the guide, ( ) Vx indicates the average over the 
measured distribution of v x and the normalization con- 
stant is chosen to satisfy YlnPn = !• 

(Ill) At the exit of the guide, the wave function is 
changed because the ceiling suddenly disappears. From 
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FIG. 2. (color). The resulting probabilities for eigenstates 
p n (a) and the corresponding Wigner function W(z,v z ), the 
sum of the Wigner functions n < 50 with the weights of the 
probabilities for the eigenstates (b). W(z,v z ) for v~ > is 
shown here because W(z, v z ) — W(z, —v z ). The values in the 
phase space point in (b) denoted by colors are in arbitrary 
units. In both figures, the best fit parameters are used. 



the sudden approximation, the wave function is contin- 
uous at the time of the sudden change as J2 m ^m^m = 
2^2 n a "ni'n where the left (right) hand side corresponds to 
the wave function with (without) the ceiling. Using the 
properties of a complete orthonormal system, we obtain 

where the probabili- 

m | . In this equation, 



/„ dz ipniztymiz) 
= \a n \ 2 and p m = 



Pn / , m Pro 

ties satisfy p n 

the cross terms are cancelled due to random phases. The 
resulting probabilities of the eigenstates p n are shown in 
Figure [2 (a). 

(IV) The detected position on the detector after the 
magnification depends not only on the height z but also 
on the vertical velocity v z at the end of the guide. Hence 
we use a kind of probability density in phase space given 
by the Wigner function [11[. The Wigner function is 
defined as 

W(z, Pz ) = ^ |°°d^*(2-^)^+K)exp(-^ 

where z is the height, p z is the momentum along z and ip 
is the wave function. The Wigner function in z-v z phase 
space can be obtained by replacing p z with mv z . The 
Wigner function has properties of J_ oc dp z W(z,p z ) = 
\tp(z)\ 2 and f^dz W(z,p z ) = \<p(Pz)\ 2 where ip(jp z ) is 
the momentum space wave function. The Wigner func- 
tion of a mixed state p = l^ jn Pn\'4'n){'4'n\ is a simple sum 
W{z,p z ) = z~2nPnW n (z,p z ) where W„ is the Wigner 
function for the n-th eigenstate. The Wigner function 
obtained for this system is shown in Figure [2] (b). We ap- 
proximate the correspondence of the phase space point at 
the guide end and the detection point on the detector by 
the classical trajectory. This approximation is supported 
by the fact that if the potential of the system has terms 
only up to the second order in the position, the motion of 
the Wigner function in the phase space can be obtained 
by the equation of motion in classical mechanics [l2[ . In 
this case, the potential has only first order term of mgz. 




TABLE I. Measured systematic uncertainties (with la). 
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FIG. 3. (color). The distributions of UCNs at the pixelated 
detector for the data (points with errors) and for the predic- 
tion with the best fit parameters (lines) for the whole region 
(a) and the lower ^'region (b). n < 50 eigenstates are con- 
sidered in the calculation. The distributions with the selected 
eigenstates are also shown in (b). 



The whole phase space was divided into a mesh with a 
size of Az x Av z =0.1 (im x 0.1 mm/s and each mesh 
point is weighted by the Wigncr function W(z,v z ). 

The predicted distribution was fitted to the data using 
a binned maximum likelihood method. Six parameters 
were used in the fit, 9, Zq, d 7 7, (3 and s. 9 denotes the 
rotation of the pixelated detector in the detector plane. 
The position of the data Z is rotated by 9 to be Z — >• Z' . 
Zq is the offset for the position of predicted events .Zprod- 
The relation of the coordinates is Z' = Z wct \ + Zq. d is 
the difference between the actual and design heights of 
the pixelated detector, hence d effectively modifies the 
magnification power, d > denotes that the actual posi- 
tion is higher than the design, corresponding to a higher 
magnification. 7 and /3 are the parameters describing 
losses inside the guide, s is the signal fraction. A flat 
background in Z is assumed. The predicted distribution 
is normalized to have the same sum of weights as number 
of events in the data. 

Events within the whole expected region 0.4 mm < 
Z' < 3.6 mm are selected and used to fill a histogram 
with 400 bins of width 8 /im. The eigenstates with n < 50 
are used in the calculation. The best fit parameters are 



s = 0.95 



+0.01 

-0.02' 



Z = -1.006 ± 0.002 mm, d = -0.03 ± 
0.01 mm, 9 = -1.2540ig;ggg| deg, 7 = 5.0±0.9x 10 4 s" 1 
and j3 = 0.07 ± 0.03. The x 2 over degrees of freedom 
(NDF) x 2 /NDF = 402.6/394. The distributions of the 
data and the prediction of the best fit are shown in Figure 
[3) The small Z' region of the distributions are shown in 
Figure [3] (b). The modulation of the prediction is mostly 
due to eigenstates with n < 15. The form of the several 
modulations observed in the data is in good agreement 
with the best fit prediction. 



Source 


AZ 


Az 


v x dispersion 


1.2 /mi 


0.1 /im 


Aberration of rod 


4.8 /xm 


0.3 /im 


Spatial resolution of detector 


3.4 jam 


0.2 /mi 


Roughness of detector surface 


10.6 /im 


0.6 fira 


Total 


12.1 /im 


0.7 /im 



The estimated systematic uncertainties are summa- 
rized in Table HI AZ is the uncertainty of Z and Az 
that of z. These uncertainties are calculated for neutrons 
with z = and v z = where the magnification power is 
the minimum value 16.5. The total uncertainty of z is 
Az = 0.7 /im. No detector effect was found which could 
fake the observed modulation. The predicted modulation 
is also robust against varying the values of fit parameters. 

In conclusion, we observed the spatial distribution of 
gravitationally bound states of UCNs using a dedicated 
novel technique. The vertical distribution of UCNs was 
magnified by a cylindrical rod and detected by a pixe- 
lated detector. The measured UCN distribution on the 
pixelated detector can be derived using the Wigner func- 
tion. The shape of several peaks of the modulation of the 
UCN distribution in the lower Z' region is in good agree- 
ment with the prediction. This is the first observation of 
gravitationally bound states with submicron resolution. 
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